A Nearly Scale Invariant Spectrum of Gravitational Radiation from Global Phase 

Transitions 



Katherine Jones-Smith^, Lawrence M. Krauss^'^, and Harsh Mathur-'^ 
^ CERCA, Department of Physics, Case Western Reserve University, 

10900 Euclid Avenue, Cleveland OH 44106-7079 
^ Also Department of Astronomy, Case Western Reserve University 

Using a large N sigma model approximation we explicitly calculate the power spectrum of grav- 
itational waves arising from a global phase transition in the early universe and we confirm that it 
is scale invariant, implying an observation of such a spectrum may not be a unique feature of infla- 
tion. Moreover, the predicted amplitude can be over 3 orders of magnitude larger than the naive 
dimensional estimate, implying that even a transition that occurs after inflation may dominate in 
Cosmic Microwave Background polarization or other gravity wave signals. 



In 1992 yj it was proposed on dimensional grounds 
that a scalar field which undergoes a symmetry break- 
ing phase transition can give rise to a scale-invariant 
spectrum of gravitational waves by virtue of causality: 
the correlation length of the field is of order the horizon 
size, as the horizon grows, uncorrelated regions come 
into contact, and the field releases energy as it relaxes. 

Following the detection of both temperature and po- 
larization fluctuations in the Cosmic Microwave Back- 
ground (CMB), there is now great interest in the pos- 
sibility of using CMB measurements to probe for a 
nearly scale invariant gravitational wave spectrum that 
is widely considered to be the "smoking gun" of in- 
flation [5]. Long wavelength primordial gravitational 
radiation would affect the polarization of the CMB |3] . 

Given the importance therefore of a possible obser- 
vation of a gravitational wave signature in the CMB, 
and proposed missions capable of observing CMB po- 
larization such as CMBPol and Planck it is imperative 
to go beyond dimensional arguments and explore the 
detailed nature of the spectrum of radiation produced 
by another mechanism such as that in [1^ from which 
the inflationary spectrum would need to be discrimi- 
nated. This is the purpose of the present work. 

Our analysis proceeds in two stages. First we ana- 
lyze the ordering of a symmetry breaking scalar fleld 
evolving in a flat FRW background. Then we compute 
the gravitational radiation produced by the relaxation 
of the scalar field. The estimates given in [T] assumed 
a two component scalar field in a Mexican hat poten- 
tial. However the two component case presents special 
complications due to the possibility that strings will 
form and frustrate the ordering. For this reason we 
work with an A^-component field and work in the sigma 
model limit. In other words we impose the constraint 
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uum expectation value of the field. Due to the con- 
straint the equations of motion for this system are 
non-linear and in general can only be solved numer- 
ically. However it has long been known in the con- 
densed matter literature on ordering kinetics that the 
problem becomes tractable in the large N limit; that 
large N method was adapted to the study of global 
phase transitions in cosmology by Turok and Spergel 
[4J whose analysis we will follow. 

A key approximation in the large N solution is to 
replace the trace of the stress-tensor with its spatial 
average in the sigma model equation of motion. Thus 
we set about solving 
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2 da d(f>a 
a dr dr 



-S/^cf,^ = -\T{T)cl>^ (2) 



where t is conformal time, a is the scale factor, and 
T(r) is the spatial average of trace of the sigma field 
stress tensor, a^T^ ^{v,t) [6j . For simplicity we first 
consider the case that the scale factor varies as a cx t^. 
The solution is constructed self-consistently. We posit 
that r -Toti'^/t^ and find 



0„(k,r) = A„r(i/2)-'5j,(fcr) 
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on N otherwise free massless fields. Here rj is the vac- 



where 0Q(k, r) is the spatial Fourier transform of the 
sigma model field, is to be determined by the ini- 
tial condition and the order of the Bessel function is 
given by 1^^ = Tq + {P - 1/2)2; ^^^^i Tq and ly will be 
determined self-consistently below. 

To model the initial symmetric state we take the 
sigma field to be completely disordered at time tq. To 
regulate the initial field distribution we introduce a 
real space cubic lattice with lattice constant d, and 
choose "white noise" initial conditions where the field 
components at each lattice site are taken to be inde- 
pendent Gaussian random variables. The variance of 
the distribution is chosen so that the constraint eq (jlj 
at each site is satisfied on average; the deviations from 
the mean vanish in the large N limit. Since the Fourier 
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amplitudes of the sigma field evolve independently ac- 
cording to eq ([3| they remain independent Gaussian 
random variables. Their correlators at later times can 
be evaluated in terms of their correlators at the initial 
time tq. 

Finally we impose self-consistency by requiring the 
average stress tensor trace T(t) computed from our so- 
lution (/)Q(k, r) to equal the postulated form — ToTy^/r^. 
This requirement determines that i' — 1 + (3 and 

In summary, in the large N limit, the two-point 
correlators of the sigma model field are given by 
(0Q(k,ri)0^(p,T2)) = {2nf6ajS{k+p)C{k,Ti,T2) and 



C(fc,ri,T2) = 8: 
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Four-point and higher correlators may be computed by 
use of Wick's theorem. 

Experience with ordering dynamics suggests for ex- 
ample that the equal-time two-point correlator should 
have the scaling form C(A:,t, r) = T"/(fcT^). Indeed 
eq Q has this form with dynamical exponent z — \ 
and a — i. That the weight of a mode, |(/)Q(k,T)p, 
grows as r'^ before it comes into the horizon (r ^ 
can be understood as a consequence of the constraint 
eq ([TJ. Initially all modes have the same weight, but 
as the universe expands, the field orders on the horizon 
scale, and short wavelength modes that have already 
come into the horizon decay. To satisfy the constraint, 
spectral weight shifts into the long wavelength modes 
that lie inside a sphere of radius 1 /r in wave-vector 
space. Since the total weight is conserved, the weight 
of these modes grows as r^. 

So far we have assumed that the scale factor varies as 
a power law a oc . This is valid deep in the radiation 
and matter dominated eras with P — I and 2 respec- 
tively. (Note that these are the relevant eras for CMB 
polarization studies.) To model the crossover between 
the two eras we take the scale factor to be 
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where r* is the crossover time. We define a time depen- 
dent exponent 6(t) = {t /a)da/dT. For a pure power 
law scale factor this coincides with the exponent /3, and 
for the crossover form, 6(r) has the asymptotic values 
1 for r <C T* and 2 for r » r* . Eq ([2]) is no longer sol- 
uble in closed form when a is given by eq ([5]) , but since 
6(t) varies slowly, it is an excellent approximation to 
work with the solution for fixed /3 derived above and 



simply make judicious replacements of ^(t) in place of 
(3. Turok and Spergel have verified the accuracy of 
this "adiabatic" approximation numerically. Note that 
in this case the sigma field correlators depart weakly 
from a rigorous scaling form since h varies slowly with 
r. 

We now turn to the calculation of the gravitational 
radiation. We decompose the metric into g^t^ = 
VfiviT) + hf^iy where 77^^ = a'^{T){dT^ - dx^ - dy"^ ~ dz^) 
is the fiat space FRW metric and h^^ is the pertur- 
bation due to the evolving scalar field. As customary 
in linearized relativity we decompose the perturbation 
into scalar, vector and traceless tensor parts via /iqo = 
2a?^, hoi — —a'^Wi and hij — —2a^Sij + 2a^tl)6ij. It is 
most convenient to use the transverse gauge diWi = 
and diSij = and to work with the spatial Fourier 
transform of the perturbation. In this representation, 
for each wave-vector there are two transverse degrees 
of freedom whose amplitudes are conventionally writ- 
ten /i+(p, r) and hxip,T). These amplitudes obey the 
linearized Einstein equations 
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and a similar equation for hx- For the wavevector 
p along the z-axis the non-vanishing components of 
the metric perturbation are related to the amplitudes 
via sii(p,r) = -S22(p,t) = h+ and si2(p,t) = 
■S2i(Pj''") = ^x- For other orientations of the wave- 
vector, Sij can be computed by rotating the compo- 
nents of the tensor appropriately. The gravitational 
modes and hx are driven by the stress in the sigma 
model field. For p along the z-axis the driving terms 
are related to the components of the sigma-field stress 
tensor via T+(p,t) = 47rG[rii(p, r) - T22(p,t)] and 
Tx(p,r) = 87rGTi2(p, t)[3. For other orientations of 
the wave- vector the source terms can be constructed by 
appropriately rotating the components of the sigma- 
field stress-tensor. 

Before we embark on a full solution it is useful to 
gain some qualitative understanding of eq In 
Minkowski space a = 1 and gravitational mode am- 
plitudes respond as simple harmonic oscillators. In the 
FRW case there is an additional time-dependent damp- 
ing produced by the expansion of the Universe. For 
a (X and in the absence of external driving (r+ = 0) 
it is easy to show that the oscillations of the gravi- 
tational wave amplitude decay as 1/a for late times 
r ^ l/p; in other words, modes redshift and decay 
once they come into the horizon. For the case that a 
mode is driven by a power law source, cx r", it is 
easy to show that for short times, r ^ 1 /p, a previously 
undisturbed gravitational mode grows as cx r"'''^. 
This is because for short times the mode amplitude be- 
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haves essentially like a driven free particle. For the case 
of interest, where the source is the sigma field stress, 
we expect that r+(p,r) « in the long time limit 
T ^ 1/p since the sigma field is ordered on sub-horizon 
scales. In other words a gravitational mode should 
evolve freely, without external driving, once it comes 
into the horizon. For short times, t <C 1/p, the driving 
term T_|_(p,t) l/v^! this follows simply from the 
scaling form of the sigma field correlators and power- 



the product pr (as can be seen by suitable rescaling). 
Fig 1 shows the strain spectrum evaluated in this way 
for the case of matter domination (3 = 2. For small pr 
we have derived the asymptotic result 
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counting. Thus we infer that |/i+(p, r' 



P cx T'^ for short 



times, T ^ l/p. 
Finally for a oc 



For long times, |ft-+(p, r)| c>c 1/a 
it is easy to see that, since the 



sigma field correlators have a scaling form, so too will 
|ft.+ (p,T)p. This translates into the observation that 
the gravitational strain spectrum (defined below) is not 
a separate function of p and r but of the product pr. 
Thus the spectrum will be scale invariant in the sense 
the term is commonly used in the literature. 

The linearized Einstein eq ([6| is soluble by conven- 
tional Green's function methods jSj: 
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For a (X the Green's function is given by 
Gp (r, ri ) = pTTi [jp^ 1 (pri ) i (pr) -n^_ i (pn ) 



here jp-i is the spherical Bessel function and 
spherical Neuman. For a given by eq ^ we evalu- 
ate the Green's functions using the adiabatic approxi- 
mation, replacing (3 with 6(t) and 6(ti) appropriately. 
The gravitational strain spectrum is defined as 
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consistent with the behavior deduced on physical 
grounds above. The constant Cfj, which gives the ap- 
proximate power at pr = 1, is w 2500 for /? = 1 and 
1000 for (3 = 2. For large pr the numerically calculated 
fall off is consistent with the expected form but 
we have not yet confirmed this by asymptotic analysis. 
Right at the peak (which occurs at pr « 3.7) we find 
the power goes as nG'^r/'^/N, in agreement with the 
estimate of ref [Ij, but factors of 27t remarkably con- 
spire to make the numerical factor k = 11,600. Thus 
the gravitational radiation produced by this mecha- 
nism is over four orders of magnitude larger than was 
estimated there. Modes that enter the horizon in the 
matter dominated era satisfy pr* <ti 1; for these modes 
Fig 1 is a very good approximation for any t ^ t* , i.e., 
the entire relevant history of the Universe. 

Our key findings are summarized in Figs 1 and 2. Fig 
1 shows the power spectrum of gravitational radiation 
P{p, t) calculated in an approximation that is valid in 
the matter dominated era. Here p is the magnitude of 
the gravitational wavevector and r is conformal time. 
As we have shown, the matter dominated era the power 
is a function of the product pr. Fig 1 shows that for 
modes of fixed p the power peaks at r '--^ l/p when 
the modes come into the horizon. For later times the 
modes redshift and consequently the power falls off as 
^-'h/r*. For earlier times the power grows as r^. Fig 1 
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Eqs Q, Q and ([s]) provide all the ingredients re- 
quired to compute the gravitational strain spectrum. 
Performing this calculation, making use of Wick's the- 
orem to evaluate the required four-point sigma field 
correlator, we arrive at our final result 
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This result simplifies if a oc t^: then b —> P and we 
may set r = 1 since the power becomes a function of 



may also be interpreted as a plot of power vs. wave 
vector for fixed time. Then, the power is seen to peak 
at wavevectors on the horizon scale p = 1/t and to 
drop off larger and smaller scales. 

Fig 2 shows a plot of the exact power, eq (10 1, as a 
function of conformal time r at wave- vector p = 1/r, 
over a much more complete range of time. Fig 2 shows 
that the spectrum is indeed scale invariant over the 
range of times r 3> t* and that it approaches the value 
expected for pure matter domination, with a jump at 
matter-radiation crossover time, and a plateau for t 
T* at the value for pure radiation domination. 

These results not only confirm expectations based on 
dimensional grounds, they provide significant motiva- 
tion to further explore in detail the possible observa- 
tional consequences of the gravitational wave spectrum 
predicted here. The surprisingly large magnitude of the 
power in gravitational waves produced by this mech- 
anism does not require inflation not to have occurred 
for the spectrum generated by such a global transi- 
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FIG. 1: Plot of gravitational radiation power vs pr where 
p is the wave- vector and r is the conformal time, computed 
assuming the scale factor appropriate for matter domina- 
tion. The plot may be interpreted to show the variation 
in power with time for modes of a fixed wave-vector or the 
variation in power with wave-vector for fixed time. The 
points are obtained by numerical integration of eq ( 10 ). The 



dashed curve is a guide to the eye. The solid line shows an 
asymptotic approximation Cplpr)^ valid for small pr. The 
numerical data are consistent with a 1/(pt)* tail for large 
pr. The peak value 11,600 occurs at pr « 3.7. 
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FIG. 2: Plot of the gravitational power vs r/r* at wavevec- 
tor p = 1/t. r, is the time of crossover from radiation to 
matter domination. The spectrum is flat for r ^ r* and 
again for r ^ r, ; thus except at the crossover, it is essen- 
tially scale invariant. The points are calculated by numer- 
ical integration of eq ( |10| l; the dashed curve is a guide to 
the eye. 



tion to be observable. Even a phase transition below 
the scale of inflation might produce a signature that 
could compete with or overwhelm that generated dur- 
ing inflation. As a result some observable difference in 
the resulting polarization signal will be necessary if we 
hope to claim such a signal may provide unambiguous 
evidence that inflation occurred in the early universe. 
Work is now in progress to calculate the CMB polariza- 
tion signature due to this radiation in order to seek out 
possible distinctions in the signatures for the different 
alternatives. 

Another figure of merit useful for experimental pur- 
poses is Ogw(p, Tnow), the present energy density per 
log frequency interval of subhorizon modes normalized 
to the critical density. It is given by ^gwiP,Tnovr) = 



p^P(p, Tnow)/6i?^ where H is the present Hubble con- 
stant. At frequencies that can be detected by ad- 
vanced LIGO, LISA and millisecond pulsars, both in- 
flation and this mechanism predict a flat spectrum for 
^gw(P7 Tnow), with relative strengths given by: 
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where — Vq is the energy density during inflation. 
Once again, depending upon relative energy scales, this 
signature could dominate over that due to inflation. 

As a final note, while we expect that the general fea- 
tures of the resulting spectrum are not dependent on 
the large N approximation which allowed an exact an- 
alytical form to be determined, it may also be desirable 
to study the case N — 2 and = 3 cases. Here the for- 
mation of topological defects (strings and monopoles) 
will also infiuence the dynamics of the symmetry break- 
ing field, leaving an imprint on the spectrum of gravi- 
tational radiation. 

We acknowledge discussions with Irit Maor and Tan- 
may Vachaspati. 
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